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The nonlocal version of the SU(2) x SU(2) symmetric four-quark interaction of the NJL type 
is considered. Each of the quark lines contains the form factors. These form factors remove the 
ultraviolet divergences in quark loops. The additional condition on quark mass function m(p) ensures 
the absence of the poles in the quark propagator(quark confinement). The constituent quark mass 
m(0) is expressed thought the cut-off parameter A, m(0) = A = 340 MeV in the chiral limit. These 
parameters are fixed by the experimental value of the weak pion decay and allow us to describe 
the mass of the light scalar meson, strong decay p — > nn and D/S ratio in the decay ai — » pn in 
satisfactory agreement with experimental data. 
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(NJ . I. INTRODUCTION 

wo: 

The effective meson Lagrangians obtained on the basis of the local four-quark interaction of the Nambu-Jona- 
<^ Lasinio(NJL) type satisfactorily describe low-energy meson physics 0, 0, y, B- However, these models contain 
ultraviolet (UV) divergences and do not describe quark confinement. Therefore, additional regularization is necessary 
in these models. Besides, it is impossible to take into account the dependence of the amplitude of different processes 
. on large external momentum in order to provide quark confinement. This situation restricts the predictive power 
of these models. Satisfactory results in these models can be obtained only for light mesons and interactions at low 
£SJ , energies in the range of 1 GeV. In order to overcome these restrictions, it is necessary to consider nonlocal versions 
t"***" ■ of these models which allow us to remove UV divergences and describe the quark confinement. 

A lot of models of this type were proposed in the last few years. Unfortunately, we cannot give here the full list of 
references concerning this activity. Therefore, we will concentrate only on the direction connected with the nonlocal 
quark interaction motivated by the instanton theories 0> El 0> • Recently, a few nonlocal models of the this type 
were proposed [^. llCtlTlLll2| . In these models the nonlocal kernel is taken in the separable form where each quark line 
contain form factor following from instanton theories. These form factors naturally remove UV divergences in quark 
loops. Thus, in |3. Q(| a nonlocal form factor was chosen in the Gaussian form f(p) = exp(— p 2 /A 2 ) where A is the 
cut-off parameter [2 1| - In |lllll2j it was proposed to use an additional condition for the form factor f(p) (quark mass 
CIh! function m(p), respectively) which lead to the absence of the poles in the quark propagator. Namely, it is supposed 
1} ■ that the scalar part of the quark propagator is expressed through the entire function 

2 7 (j t exp(VM 2 ), (1.1) 

p z + m z (p) fi 

where p is an additional arbitrary parameter. The similar condition providing confinement was used in |l3l Il4 . ITfll |. 

In this work an analogous condition will be used . However, we will take into account that each quark line contain 
square of the form factor which is expressed through the quark mass 

m(p)f 2 (p) m 2 (p) . 2 2 . 

T~ — 27T ^ TT — 27T = ex P \~P / A ) • ( L2 ) 

pA _|_ m Z pA _|_ m Z ^pj 

As a result, we obtain a simpler solution for the mass function than in [TTIIT^. In our model m(0) and the cut-off 
parameter A have a simple connection in the chiral limit m(0) = A; the function m(p) contains only one arbitrary 
parameter. We fix this parameter by weak pion decay. Then, for = 93 MeV we have m(0) = A = 340 MeV. This 
leads to reasonable predictions for the scalar meson mass, width of the decay p — > 7T7t and D/S ratio in the decay 
a% — > pir, where D,S are the partial waves of this decay . 

The paper is organized as follows. In Sect^J we consider a nonlocal four-quark interaction and after bosonization 
derive the gap equation for dynamical quark mass. The additional condition for this mass allows us to provide the 
quark confinement. In Sect lllll the masses and couplings of the scalar and pseudoscalar mesons are obtained and 
the main parameters of the model are fixed. The decay width a — > irir is calculated. In Sect II VI the vector and 
axial-vector sectors of the model are considered. The a\ meson mass, decay widths p — > tttt, a\ — > pir are calculated. 
D/S ratio in the decay a\ — > pix is estimated. The 7r — a\ mixing is studied. The discussion of the obtained results 
and comparison with others models is given in the last section. 
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II. SU(2) x SU(2) NONLOCAL QUARK INTERACTION 

The SU(2) x SU(2) symmetric action with the nonlocal four-quark interaction has the form 



S(q,q) = d 4 x {q(x)(id x - m c )q(x) 



+^ {J a Ax)j*{x) + J a {x)Ux)) + jlt(x)j^(x))Y (2.1) 

where q{x) = (u(x),d(x)) are the u and d quark fields, m c is the diagonal matrix of the current quark masses, G\ is 
the coupling constant of the scalar and pseudoscalar quark currents, G 2 is the coupling constant of the vector and 
axial- vector quark currents. The nonlocal quark currents Ji(x) are expressed as 



Ji(x) = J d x x d x 2 f(xi)f(x 2 ) q(x - Xi) Tj q(x + x 2 ), 
where f(x) are the nonlocal functions. In l|2.2|l the matrices T/ are defined as 

v a = i,rj = i7 V\r£ Q = yr-.r^ = 7 V^ Q , 

where r° are the Pauli matrices and 7 M ,7 5 are the Dirac matrices. 

In this article, we mainly consider the strong interactions. The electroweak fields may be introduced by gauging 
the quark field by the Schwinger phase factors (see, c.f. |sl. lij Ho| 

After bosonization the action becomes 

S(q,q,a,w,p,a) = J d A x {q(x)(id x - m c )q(x) - ±- (ir a {x) 2 + a{xf) + ((^») 2 + K») 2 ) 

+ J d 4 x 1 d i x 2 f(x - x 1 )f(x 2 - a:)g(as x ) (■K a (x)i 1 \ a + a(x) + p» a {x)^r a + OhV^) q(x 2 )} , (2.2) 

where a, n, p, a are the er, n, p, a\ meson fields, respectively. The field a has a nonzero vacuum expectation value 
< a >o— Co ^ 0. In order to obtain a physical scalar field with zero vacuum expectation value, it is necessary to shift 
the scalar field as a = a + <jq. This leads to the appearance of the quark mass function m(p) instead of the current 
quark mass m c 

m(p) =m c + rridynip), (2.3) 
where nid y n{p) = — &of 2 (p) is the dynamical quark mass. From the action, eq. 1)2. 2[1 . by using 

/ 

one can obtain the gap equation for the dynamical quark mass 



m dyn ( P ) = G 1 ^( P ) J d%knk) WT ^ m . (2.4) 

The right-hand side of this equation is the tadpole of the quark propagator taken in the Euclidean domain. Equations 
(|2.3|) . I|2.4I) have the following solution: 

m(p) = m c + (m q - m c )f 2 (p), (2.5) 

where m q = m(0). 

In order to provide quark confinement we propose the following anzatz for the quark mass function m(p). We 
suppose that mass satisfies the following condition in the chiral limit 



m (p) 



m? (p) + p 2 



= exp (-p 2 /A 2 ) . (2.6) 
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The form of the left-hand side of this equation coincides with the integrand in the gap equation i|2.4[l . From eq. (|2.6() 
we obtain the following solution [22|: 

m(p)=( — jjL, X /2 \ (2-7) 



exp (p 2 /A 2 ) - 1 

here we have only one free parameter A; m(p) does not have any singularities in the whole real axis and exponentially 
drops as p 2 — > oo in the Euclidean domain. From ea. H2.5jl follows that the form factors have a similar behavior that 
provides the absence of UV divergences in our model. At p 2 = the mass function is equal to the cut-off parameter A, 
m(0) = A. The pole part of the quark propagator also does not contain singularities that provide quark confinement 

1 1 - cxp (- P 2 /A 2 ) 



m 2 (p) + p 2 p 2 
When taking into account the current quark mass ea. l|2.6f) can be modified as follows 

m 2 (p) — m 2 



m 2 (p) + p 2 



exp(-(/ + m 2 )/A 2 ). (2.9) 



Here m 2 is introduced in the form that conserves the analytical properties of the mass function m(p). Then the mass 
function takes the form 

/ m 2 c +p 2 e Xp (-(p 2 + m 2 c )/A 2 ) \ 1/2 

m( P> { 1 -exp(-(p2 + m 2)/ A 2) J [Z - W) 



III. PSEUDOSCALAR AND SCALAR MESONS 



Let us consider the scalar and pseudoscalar mesons. The meson propagators are given by 

where Mc^ are the meson masses, g<7,Tr(p 2 ) are the functions describing renormalization of the meson fields and 
n CT ,7r(p 2 ) are the polarization operators defined by 



n ^ (p2) = l (|^ J d 4 fc/ 2 (fc 2 )/ 2 (4)s P [S(fc_)r^s(fc + )r ff ^] 



(3.2) 



where k± = k ± p/2. 

For calculation of these integrals it is necessary to rewrite these expressions in the Euclidean space where the 
form- factors (and quark masses) are the exponentially decreasing functions. Then ea. H3.2jl takes the form: 

a 2A^ e f 4 P a ^{k\p 2 ,p-k) 
a ' w[P ' ~ (2^) 4 m2 J aEK (k 2 + +m(k 2 + ) 2 )(k 2 _+m(k 2 _) 2 ) { ' 

The functions P ai7r (k 2 ,p 2 ,p-k) are the Dirac trace multiplied by m(fc+),m(fc_). In ea. (|3.3|l all momenta are Euclidean. 
In the description of the meson properties it is necessary to make the analytical continuation of this expression over 
external momenta p to the Minkowski space. Let us emphasize that at our anzatz for a quark mass function only the 
functions P CTj7r (fc 2 ,p 2 ,p- k) contains nonanalytical root cuts, whereas there is no problems with analytical continuation 
of the denominator. 

The meson masses are found from the position of the pole in the meson propagator 

U^Ml J = Gl\ (3.4) 

and the constants g a ^{M 2 ^ are given by 

<a«J = ^^U=^. (3.5) 
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Firstly, let us consider this model in the chiral limit. The pion constant 3^(0) is not depend on parameter A and 
takes the form 



here £ is the Riemann zeta function. 
The gap equation has the simple form 



, 2tt 2 

= — . (3.7) 



The quark condensate is 

oo 

(m)o = -^[duu-^ n . (3.8) 
47r z J u + m [u) 
o 

The Goldberger-Treiman relation is fulfilled in the model of this kind 0, 0, 0, 0] 



m,. 



Fa = (3.9) 

From eas.(pO) )l .lpn? )l the value A = m q = 340 MeV is obtained for = 93 MeV. Then, from eos.lpTTjl. (|3~%|) we obtain 

Gi = 56.6 GeV, (qq) = -(188 MeV) 3 . (3.10) 

In the description of pion mass it is necessary to introduce the nonzero current quark mass m c . In our model 
M 2 <C A 2 . Therefore, we can consider only the lowest order of the expansion in small p 2 . Then, one gets from ea. (|3.1|) 

oo 



Ml=gl(Q)\G^-^- 2 / dun J ^> . (3.11) 

By using the expression for G\ from the gap equation l|2.4[l . the Gell-Mann-Oakes-Renner relation can be reproduced 

Ml = -2^0^ + O{ml). (3.12) 

From ea. l|3.12|l with M v = 140 MeV we can estimate the value of the current quark mass m c w 13 MeV. Other model 
parameters in this case change very little 

A = 343MeV, g w (M w ) — 3.57, Gi = 56.5GeV, (qq) = -(189 MeV) 3 . (3.13) 

Therefore, in calculations of the amplitudes of various processes we can use the values of parameters taken in the 
chiral limit. 

With the help of the parameters l|3.10[) we get for sigma meson M a — 420 MeV and g a (M a ) — 3.85. The amplitude 
of the decay a — > nn is equal to A^^+tt-) — 1-67 GeV. Then, the total decay width is r^^^ = 150 MeV. 
Comparing these results with experimental data one finds that M a is in satisfactory agreement with experiment 
M ex P _ 400 _ 120 o ; however, the decay width is very small Tf*> = 600 - 1000. 

IV. VECTOR AND AXIAL- VECTOR MESONS 

The propagators of the vector and axial-vector mesons have the transversal and longitudinal parts 

D HV rpiiUjjT , L ^ D L u I) 

where = g» y - p^p u jp 2 , = p^p v jp 2 and 



D T = - = 9 P' a ^ D L = - ( 4 ■> i 

p ' ai g^ + uj^) m^-p 2 ' "> ai G^ + n^y 
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Here, nj ai and ITf ai are the transversal and longitudinal parts of the polarization operator ^^ v ai {p 2 ) 

on r 

K:aM) = l j^i J d 4 fc/ 2 (fc_)/ 2 (fc+)s P [s(fc_)r p , ai s(A :+ )r,, ai ]. 

The constant G2 is fixed by the p-meson mass 

g 2 - x = -n T p (M p ) 

and G2 = 6.5 GeV -2 . Then the ai-meson mass is equal to 970 MeV. 
The constants g p . ai (M 2 0i ) are equal to 

dU T (v 2 ) 

9-^Ml, ai ) = -^^\ P ^M lai . (4.3) 

From eq.gSJ we obtain g p {M p ) = 1.23, g a (M ai ) = 1.33. At p 2 = we have g p (0) w 2, g o (0) « 2.5. (see also FigHJ). 
The decay p — > 7T7t is described by the triangle quark diagram. The amplitude for the process is 

^ P ^=a( P ^)(Qi-Q2r (4-4) 

where qi are momenta of the pions. We obtain a( p ^ 7r7r ) = 5.72 and the decay width Y( p ^ vv \ = 135 MeV which is in 
qualitative agreement with the experimental value 149.2 ± 0.7 MeV (l7j| . 

The decay a\ — > pn is described in a similar manner. The amplitude for the process a\ — > pn is 

A Ta^) = Ha^)9^ + h*^)P V ^ (4.5) 

where p, q are momenta of ai,p mesons, respectively. We obtain a( Ql ^ p7r ) = 2.68 GeV, b( ai ^ p ^ = 16.71 GeV -1 . 
Amplitude of the decay ax — > pn contains D and S waves. The ratio of these waves has the form(see [Tol fl5| ): 

/S ~ ^(E p + 2M p )a {ai ^ )+ b ia ^ p7r) M ai \t\ 2 - °-° 6 ' (46) 
|t I 2 = A(Af 2 i , M 2 , M 2 )/(2M ai ) 2 , E 2 = M 2 + |^| 2 , A(a, b, c) = a 2 + b 2 + c 2 - 2ab - 2ac - 2bc 

This ratio is in satisfactory agreement with experimental data D/S cxp — —0.108 ± 0.016. The decay width equals 
r( Q1 ^ P 7r) = 90 MeV. This value is noticeably smaller than experiment 250 — 600 MeV [l7| . 

The longitudinal component of the ai-meson field is mixed with the pion. The amplitude describing this mixing 
has the form 

where p is the momentum of the pion, and C( ff _ >ni )(0) in the chiral limit is equal to 190 MeV. The additional pion 
kinetic term from the ir — a\ mixing is ALkin — A ■ p 2 ir a (p) 2 /2. In the chiral limit A is equal to 

C 2 (0) 

A - ( "^ ai) wC? w „ ,(0)G 2 / ff 2 (0)«0.04. (4.7) 

A is small; therefore, the effect of the n — a\ mixing can be neglected. 

V. DISCUSSION AND CONCLUSION 

In this work we have considered a possibility of constructing the SU(2) x SU{2) symmetric nonlocal chiral quark 
model providing the absence of UV divergences and quark confinement. These features of the model are specified by 
the nonlocal kernel which appears in the four-quark interaction. Such a structure of the four-quark interaction can 
be motivated by the instanton interactions [|| 0, @ • 

The pseudoscalar, scalar, vector and axial-vector mesons have been considered in the framework of this model. 
The masses and strong coupling constants of the mesons were described. It was shown that the functions describing 
renormalization of the meson fields noticeably decreased at large p 2 in the physical domain(see Fig^). 
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Among satisfactory predictions of the model are the mass of the cr-meson , the decay width p — > tttt and the D/S 
ratio in the decay a\ — * pit. 

However, in the description of the a\ meson mass and decay widths a — > ttt: , a± — > pn our results are noticeably 
smaller than experimental data. Note that the width of the decay a\ — > pn strongly depends on mass of the a\ meson. 
Indeed, for M ai — 1.26 GeV we have r^ ai ^, pn -j w 200 MeV that is in qualitative agreement with experiment. 

It is useful to compare the obtained results with the analogous results obtained in the local NJL model 2] and 
other nonlocal models with quark interaction of separable type [Tol Il5| . 

Remind that in the local NJL model the cut-off parameter A^ NJL ^ = 1.2 GeV and the constituent quark mass 
m = 280 MeV are used. These parameters are fixed by the decays tt — > nu{f n — 93 MeV) and p — » Trir(g p = 6.14). 
These parameters lead to the quark condensate (qq)o — —(293 MeV) 3 and current quark mass m c — 3 MeV. In the 
present model m = 340 MeV plays the role of the constituent quark mass, whereas our parameter A = 340 MeV 
corresponds to the effective cut-off parameter A cff ss 800 MeV. As a result, we obtain (qq)o = —(188 MeV) 3 and 
m c = 13. Remind that these values correspond to the physical pion mass. 

Let us consider 7r — a\ mixing in these models. In the local NJL model the amplitude describing the -k — a\ mixing 
equals A^^f a ^ = iy/ftmp* 1 . Therefore, the coefficient C^f^K in the equals 680 MeV. This value is 3.5 times larger 
than in the present model. As a result, it leads to the noticeable additional renormalization of the pion field in the 
local NJL model gf JL) = gi NJL) ■ Z 1 / 2 = m/%, where Z = (l - 6m 2 /M^) -1 w 1.4 and gi NJL) = gf JL \ in our 
model Z = 1.04. Therefore, in the local NJL model the ir — a\ mixing plays a more important role. 

Let us compare also the amplitude of the decay width a — > -kit in these models. In the local NJL model without 
taking into account the 7r — a% mixing in the external pion legs this amplitude equals A^ N _^ L } _ = Amgi NJL ^ Z 1 / 2 = 4 
GeV. This amplitude is 2.4 times larger than in the present model. However, after taking into account the it — a\ 
mixing this amplitude takes the form A\_^ + _ = Amgi N JL ' ) z^ 3 ^ 2 = 2 GeV. This is close to the amplitude obtained 
in our work. This leads to a noticeable decrease in the decay width T^^,^ — 190 MeV which also becomes smaller 
than experimental data and is in qualitative agreement with the result of our model. The mass of the a meson in 
NJL is 570 MeV and approximately 30 % larger than the result obtained here. Both the values do not contradict the 
experimental data. 

The decay p — > nn in |2| is used for fitting model parameters while in our model we predict it. The mass of the 
a\ meson obtained in our model practically coincides with the mass predicted in the local NJL model, Ma^ JL ^ ~ 1 
GeV. 

In what follows we would like to compare our result with the nonlocal model ^(| . In this model a similar separable 
instanton-motivated form of the interaction is also used. The main difference of our model with that of [l(| is connected 
with an additional requirement on a quark propagator providing quark confinement. The quark mass function in our 
model contains only one arbitrary parameter instead of two parameters in |l0j . In spite of less freedom in choosing 
model parameters, our results is close to the results obtained in 0(see table 1). 

It is interesting also to compare our results with those obtained in |l5j . where the quark propagator is expressed 
through the entire functions which are similar to the function used in our work (see ea. l|2.6f) ). The quark interaction 
in this work has a separable type which is obtained from the quark-gluon interaction with the help of a modified gluon 
propagator. In this work, the decays p — * 7T7t and a\ — > pit also have been calculated. The decay ratio D/S is close 
to ours while the decay widths strongly differs (see table 1). 

The failure of the local NJL model and its nonlocal extensions to describe the cr-meson is expectable. Similar 
problems appeared in the QCD sum rule method. In the scalar channel with vacuum quantum numbers the corrections 
from different sources may be valuable. Indeed, it has recently been shown that the 1/N C corrections in this channel 
are rather big 18j, and the Hartree - Fock approximation may be inadequate in this case. Moreover, for a correct 
description of the scalar meson it is necessary to take into account the mixing with the four-quark state [l9| and the 
scalar glueball 20] . 

In future, we plan to describe electromagnetic interactions in the framework of this model, calculate the e.m. pion 
radius, polarizability of the pion and consider the processes 7r° — > 77, 7* — > jir in a wide domain of photon virtuality. 
We also plan to generalize this model to the U(3) x U(3) chiral group by introducing new parameters: mass of the 
strange quark m s and the cut-off A s which allows us to describe intrinsic properties and interactions of strange mesons. 

The authors thank A. E. Dorokhov for fruitful collaboration, and D. Blaschke, C. D. Roberts and V. L. Yudichev 
for useful discussions. The work is supported by RFBR Grant no. 02-02-16194 and the Heisenberg-Landau program. 
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FIG. 1: Momentum dependence of the mesons renormalization functions. 



Quantity 


Our model 


[2] 


i 


PJ 


[15] 


[17] 


M a (MeV) 


420 


570 


443.2 


465.8 






400 - 1200 


r ff _™ (MeV) 


150 


190 


108 


135.1 






600 - 1000 


T p ^ (MeV) 


135 


150 


126 


114 


356 


259 


149.2 ±0.7 


M ai (MeV) 


970 


1030 


946.8 


1061.5 


1340 


1230 ±40 


Y ai ^ (MeV) 


90 


290 


44 


376.2 


4020 


385 


250 - 600 


D/S 


-0.06 




-0.048 


-0.087 


-0.092 


-0.075 


-0.108 ±0.016 



Table 1. The comparison of the physical results are obtained in local and nonlocal quark models. In local NJL model 
the decay width p — > ttt: is used for fitting of the model parameters. Two set of values corresponding to the different 
choice of model parameters are given in column 4, 5 (see table 4 in and table I in [T^l. 



